Supernova constraints on decaying vacuum cosmology 
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There is mounting observational evidence that the expansion of our Universe is undergoing a late- 
time acceleration. Among many proposals to describe this phenomenon, the cosmological constant 
(A) seems to be the simplest and the most natural explanation. However, despite its observational 
successes, such a possibility exacerbates the well known A problem, requiring a natural explanation 
for its small, but nonzero, value. In this paper we consider a cosmological scenario driven by a 
varying cosmological term, in which the vacuum energy density decays linearly with the Hubble 
parameter, A oc H. We show that this A(t)CDM model is indistinguishable from the standard 
one (ACDM) in that the early radiation phase is followed by a long dust-dominated era, and only 
recently the varying A term becomes dominant, accelerating the cosmic expansion. In order to test 
the viability of this scenario we have used the most recent type la supernova data, i.e., the High-Z 
SN Search (HZS) Team and the Supernova Legacy Survey (SNLS) Collaboration data. In particular, 
for the SNLS sample we have found 0.27 < fi m < 0.37 and 0.68 < H < 0.72 (at 2a), which is in 
good agreement with the currently accepted estimates for these parameters. 

PACS numbers: 98.80.Es; 95.35.+d; 98.62.Sb 



I. INTRODUCTION 

Over the last years, a considerable number of high- 
quality observational data have transformed radically the 
field of cosmology. Results from distance measurements 
of Type la supernovae (SNe la) 0, H 

combined with 

Cosmic Microwave Background (CMB) observations £|, 
dynamical estimates of the clustered matter and age 
measurements of the oldest structures @, seem to indi- 
cate that the simple picture provided by the standard 
cold dark matter scenario is not enough. These obser- 
vations are usually explained by introducing a new hy- 
pothetical energy component with negative pressure, the 
so-called dark energy or quintessence (for recent reviews 
on this topic, see [3]). Besides its consequences on fun- 
damental physics, if confirmed, the existence of this dark 
component would also provide a definitive piece of in- 
formation connecting the inflationary flatness prediction 
with astronomical data. 

On the other hand, from a purely theoretical view- 
point, the existence of a dark energy is related to an old 
problem of quantum field theories and theoretical cos- 
mology, namely the role of vacuum in the cosmic evo- 
lution [8|. Arguments of covariance and symmetry lead 
to an energy-momentum tensor for the vacuum of the 
form T^ v = A<?p„, where A is a scalar function which, 
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in spatially homogeneous and isotropic space-times, may 
be, at most, a function of time only. Therefore, the 
vacuum acts as a cosmological term, that is, as a per- 
fect fluid with negative pressure given (in a comoving 
frame) by p\ = —pA = —A (we work in units where 
M P ee (SttG)- 1 / 2 = c= 1). 

Nevertheless, any tentative estimation of the vacuum 
energy density by quantum field theories in flat space- 
time leads to a divergent result, and any natural cutoff 
we may choose to impose in those calculations leads to a 
vacuum contribution at least forty orders of magnitude 
bigger than the observed limits . A possible way out 
of this trouble is to postulate some cancellation mecha- 
nism that leads to an exactly null vacuum contribution, 
the observed dark energy being due to a genuine cosmo- 
logical constant or to other fields like quintessence or any 
other else. 

However, a more careful look at the problem of vac- 
uum energy may suggest another possibility. The di- 
vergent result obtained by quantum field theories in flat 
space-time cannot, rigorously speaking, be used in the 
context of curved, expanding backgrounds. On the other 
hand, in Minkowski space-time the Einstein tensor is null, 
and, therefore, any vacuum contribution to the energy- 
momentum tensor, divergent or not, should be canceled 
by a bare cosmological constant in Einstein equations. 
Now, if we could obtain the vacuum energy density 
in the expanding background, we should subtract the 
Minkowskian contribution, obtaining a "renormalizcd" 
vacuum density which would depend on the curvature, 
being very high for early times, but decreasing as the 
universe expands. 

As is well known, the Bianchi identities G^.^ = lead, 



2 



via Einstein equations, to = 0, which is an expres- 
sion of energy-momentum conservation in the presence 
of the gravitational field ^3] . In the FLRW space-time, 
this equation leads to p + 3H(p + p) = 0, where p and 
p are the total energy density and pressure, respectively, 
while H — a/a is the Hubble parameter. By introducing 
the energy densities and pressures of vacuum and mat- 
ter/radiation, one finds 



pj + 3H(p 1 + p y ) 



-A 



(1) 



The above equation is equivalent to a continuity equation 
for matter in the presence of a source —A, meaning that 
the process of vacuum decay is concomitant to a process 
of matter production, a general feature of the vacuum 
state of any non-stationary space-time. We are proba- 
bly far away from a definite theory of quantum vacuum 
in curved backgrounds or from the correspondent micro- 
scopic description of vacuum decay. Alternatively, we can 
consider effective, theoretically or empirically motivated, 
decaying laws for the vacuum density, exploring its ef- 
fects by means of macroscopic equations like Such 
an approach has an old history in the literature [t| and 
a renewed interest in recent years |Toj . 

In this regard, a viable possibility has been proposed, 
with the vacuum density decaying as A = aH, where 
a » has the order of the cube of the energy scale of 
the QCD vacuum condensation By using the ob- 

served values of m T and H, it is straightforward to verify 
that the above law provides a value for A very close to 
the value presently "observed". Naturally, the theoret- 
ical justification for this decaying law is based on some 
phenomenological hypothesis and, as such, needs to be 
verified from a fundamental theory viewpoint. However, 
the important aspect here is that it leads to cosmologi- 
cal solutions in agreement with the standard scenario for 
the evolution of the Universe, as explicitly shown in |12| . 
The above decaying law, for instance, leads to an early 
radiation-dominated phase where the vacuum term and 
the photon production are negligible, and where the scale 
factor and temperature evolve exactly as in the standard 
Friedmann solution. This phase is followed by a matter- 
dominated decelerating era, long enough to allow struc- 
ture formation, and during which the vacuum density and 
matter production are dismissable until very recently. Fi- 
nally, the Universe switches to an accelerated expansion 
driven by the vacuum, which tends asymptotically to de 
Sitter solution. If we consider the present relative matter 
density (the only free parameter of the model, besides the 
Hubble constant) around 1/3, we also obtain the present 
age parameter -ffo^o ~ U which is in good accordance 
with current observations [(|. 

In this paper we are particularly interested in testing 
the viability of the above scenarios in light of the latest 
supernova (SNe la) data, as provided recently by Riess 
et al. and by Astier et al. Q . As is well known, these 
two SNe la samples constitute the compilation of best 
observations made so far and provide the most direct 
evidence for the observed late-time acceleration of the 



universe. We also discuss other observational quantities, 
as the deceleration parameter q, the transition redshift zt 
(at which the expansion switches from a decelerated to 
an accelerated phase) , and the total expanding age of the 
Universe. In Sec. II we revise the main expressions and 
predictions of the cosmological solution with A = aH. 
The observational quantities of the model are discussed 
in Sec. III. Sec. IV presents our SNe la analysis and a 
discussion on the observational constraints. In Sec. V we 
end the paper by summarizing our main conclusions. 



II. THE MODEL 

In the context of a FLRW space-time with null spatial 
curvature, the Einstein equations lead to the Friedmann 
equation 



p = p y + A = 3H 2 



(2) 



which, together with QJ, the equation of state for the 
matter fields [p~ t = (7 — l)p 7 ], and a decaying law for A, 
completely describes the evolution of the scale factor and 
densities. As discussed earlier, we will consider A = aH, 
so that by combining the above expressions we find 



2H + 37H 2 - ajH = 0. 



(3) 



By imposing the conditions H > and p 1 > 0, we can 
also obtain the solution for the scale factor i- e -> 



a(t) = C[cxp (o-7*/2) - 1]~ 



(4) 



where C is the first integration constant and the second 
one has been set equal to zero, in order to have a = 
at t = 0. Note that the vacuum and matter densities 
are, respectively, given by A = aH and p 1 = p — A = 
3H 2 — aH. Thus, by using Eq. J3J, it is also possible to 
rewrite them as 



"■■ = tu) 



37/2 



and 



A 



a 

¥ 



c\ 37/2* 

a) 



37/2' 



(•5) 



(G) 



A. Radiation-dominated Era 

For the radiation epoch (7 = 4/3), Eq. 10} is given as 
a(t) = C [exp (2at/3) - 1] 1/2 , (7) 
so that in the limit of small times (at <C 1), we have 

a(t) sa y/2C 2 at/3, (8) 
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which is the same time dependence as in the standard 
scenario. Eqs. ©-© reduce now to 



Pr 



a 2 C 4 a 2 C 2 



3a 4 



3a 2 



and 



A = T 



3a 2 



while, in the limit a — > 0, we have 
a 2 C A 3 



Pr = 



3a 4 At 2 ' 



and 



a 2 C 2 
3a 2 



A 



a 
2t' 



(9) 



(10) 



(11) 



(12) 



From the above expressions, it is straightfoward to see 
that for small times the expansion is completely driven 
by the relativistic matter with its energy density scaling 



as a 



B. Matter-dominated Era 

For the dust phase (7 = 1), Eq. (JIJ scales as 

o(t) = C[exp(at/2) - 1] 2/3 , (13) 

which means that for small times (compared to the 
present), it can be approximated by 



a(t) = C(at/2) 2/3 . 



(14) 



Note that the above expression has the same time depen- 
dence as in the standard scenario. This, in other words, 
amounts to saying that the varying cosmological term 
starts dominating only very recently, which guarantees a 
large enough dust-dominated era. 

For this matter-dominated epoch, Eqs. ©-© read 



a 2 C 3 a 2 C 3/2 



3a 3 



3a 3 / 2 



and 



a 2 a 2 C 3 ' 2 
3 + lte 3 7 2_ ' 



(15) 



(16) 



The first term in (|15[) gives the usual scaling of non- 
relativistic matter fields, whereas the second term is re- 
lated to the production of matter at the expenses of the 
vacuum decay. Note that in the limit of large times 
(at ^S> 1), Eq. i|13|) leads to the de-Sitter solution, i.e., 



a(t) = Cexp(cri/3). 



(17) 



Note also [from i|15|) and i|16|) ] that while p m tends to 
zero A tends to a genuine cosmological constant. 



III. OBSERVATIONAL QUANTITIES 



The Friedmann equation 
epoch can be rewritten as 



(0) for the dust-dominated 



H(z) = H a 



l-O, 



n m (i + z 



,3/2 



(18) 



where il m and H stand for the current values of the 
relative matter density and Hubble parameter. From the 
above equation, it is straightforward to show that the de- 
celeration parameter, defined as q — —aa/a 2 , now takes 
the following form 



2 S ' r - 



(1 



i-n m + n m (i 

or, at the present time (z = 0), 

3 

q(z = 0) = -O rn 



;) 3 / 2 



1, 



(19) 



(20) 



Note that for any value of f2 m < 2/3 (as indicated by 
clustering and dynamical estimates [5j), the present-day 
cosmic expansion is accelerating, which seems to be in 
full agreement with current supernova observations 0, Q 
(see Sec. IV). From Eq. lfTO|) . it is also possible to obtain 
the transition redshift zt at which the Universe switches 
from deceleration to acceleration, i.e., 



ZT 



1 



-1 2/3 



1. 



(21) 



As one may anticipate, due to the process of matter pro- 
duction resulting from the vacuum decay, the transition 
redshift zt m this model will be always higher than (but 
of the same order of) the transition redshift in a ACDM 
model for the same value of f2 m . In reality, the cosmic 
acceleration in the presence of dust matter and a cosmo- 
logical term is given by 



6- = 2A- 

a 



(22) 



Therefore, the net effect of the additional terms in l|15|l 
and (|16fl is to increase the acceleration for a given value 
of p m . 

It is also possible to verify that the transition occurs 
soon after the first term in (|15Jl is supplanted by the sum 
of the second term and A. In other words, the late-time 
acceleration starts just after the end of the dust matter 
epoch, which occurs at the redshift 



(1 + V2) 



Finally, as shown in Ref. 
eter can be expressed by 



1 

ilrr, 



1 



2/3 



1. 



(23) 



[12] ]. the present age param- 



I ln(On) 



(24) 



n m - 1 ! 

which means that for the current accepted interval for 
the matter density parameter Sl m = 0.30 ± 0.05 (2a) 0, 
one finds H t ~ 1.15 ±0.08, i.e., in accordance with 
current age parameter estimates 0. 
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FIG. 1: Hubble diagram for 157 supernovae from HZS Team [Panel (a)] and 115 supernovae from SNLS Collaboration [Panel (b)]. 
As indicated in the figure, the curves correspond to a fixed value of the Hubble parameter, h (= i/o/100Km.s _1 .Mpc _1 ) = 0.7, 
and selected values of fi m . For the sake of comparison the current standard cosmologial model, i.e., a flat ACDM scenario with 
fi m = 0.27, is also shown. 



IV. SUPERNOVA CONTRAINTS 

SNe la observations are certainly among the most re- 
markable findings of modern observational cosmology 
and provide the most direct evidence for the observed 
late-time cosmic acceleration. In this Section we test the 
viability of the decaying vacuum scenario discussed above 
through a statistical analysis involving the most recent 
SNe la data, namely, the High-Z SN Search (HZS) Team 
@ and the Supernova Legacy Survey (SNLS) Collabora- 
tion data [3j- 

A. SNe la Samples 

The total sample presented by the HZS Team con- 
sists of 186 events distributed over the redshift interval 
0.01 < z < 1.7 and constitutes the compilation of the 
best observations made so far by them and by the Super- 
nova Cosmology Project plus 16 new events observed by 
Hubble Space Telescope (HST). This total data-set was 
divided into "high-confidence" (gold) and "likely but not 
certain" (silver) subsets. Here, we will consider only the 
157 events that constitute the so-called gold sample. 

The current data from SNLS collaboration correspond 
to the first year results of its planned five years survey. 
The sample includes 71 high-z SNe la in the redshift 
range 0.2 < z < 1 and 44 low-z SNe la compiled from 
the literature but analysed in the same manner as the 
high-z sample. This data-set is arguably (due to multi- 
band, rolling search technique and careful calibration) 
the best high-z SNe la compilation to date, as indicated 



by the very tight scatter around the best fit in the Hubble 
diagram and a careful estimate of systematic uncertain- 
ties. Another important aspect to be emphasized on the 
SNLS data is that they seem to be in a better agreement 
with WMAP results than the gold sample (see, e.g., 
for a discussion). In what follows we briefly outline the 
main assumptions for our analysis (see also [l4| for some 
recent SNe la analysis). 

B. Statistical Analysis 

The predicted distance modulus for a supernova at red- 
shift z, given a set of parameters s, is 

fip(z\s) = m- M = 51ogdi+25, (25) 

where m and M are, respectively, the apparent and ab- 
solute magnitudes, the complete set of parameters is 
s = (Hq, fl m ), and di, stands for the luminosity distance 
(in units of megaparsecs) , 

dL = « 1 + z )£*E&Ty (26) 

with H(x; s) being the expression given by Eq. (JTSJl. 

We estimated the best fit to the set of parameters s by 
using a \ 2 statistics, with 

x2 = £[/4(±)_aM! (27) 

where N = 157 and 115 for gold and SNLS samples, 
respectively, iA(z\s) is given by Eq. n % (z) is the 
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FIG. 2: The results of our statistical analyses. (Panel a) Confidence regions in the f2 m — h plane for the gold sample of 157 
SNe la. The best-fit parameters for this analysis correspond to fi m = 0.39 and h — 0.64. (Panel b) The same as in Panel (a) 
for the SNLS sample of 115 SNe la. In this case, the best-fit model happens at f2 m = 0.32 and h = 0.70, with \r — 1-0. 



extinction corrected distance modulus for a given SNe la 
at Zi, and cr, is the uncertainty in the individual distance 
moduli. 

Figures (la) and (lb) display the Hubble diagram 
for a fixed value of the Hubble parameter, h (= 
ifo/100Km.s _1 .Mpc _1 ) = 0.7, and some selected values 
of f2 m . For the sake of comparison our current standard 
model, i.e., a flat ACDM scenario with Q, m — 0.27, is 
also shown. Note that for the interval of Sl m considered 
the predicted magnitude-redshift relation is very similar 
in both classes of scenarios, i.e., ACDM and A(i)CDM. 
In Figs. 2a and 2b we show the results of our statistical 
analysis. Confidence regions (68.3%, 95.4% and 99.7%) 
are shown in the Q m — h plane by considering the gold 
and SNLS samples respectively. The best-fit parameters 
for the gold analysis are fi m = 0.39 and h = 0.64, with 
the reduced x 2 = Xmin/ 1 ' — 1-15 (v is defined as degrees 
of freedom). At 95% c.l. we also obtain the intervals 
0.34 < fl m < 0.44 and 0.62 < h < 0.66. Note that 
these estimated values for the matter density parameter 
are considerably large, being only marginally compatible 
with the current accepted range, i.e., Q m = 0.30 ± 0.05 
(2a) We also note that the above reduced value of 
X 2 is very similar to the one found for the flat ACDM 
scenario and is slightly larger than the one obtained for 
the ACDM model with arbitrary curvature (x 2 — 1.12). 
For the SNLS analysis (Fig. 2b), the best-fit parameters 
are fi m = 0.32 and h — 0.7 (x 2 — 1-0), which correspond 
to an accelerating universe with qo ~ —0.52, a total ex- 
panding age to — 15.7 Gyr, and a transition redshift 
z T ^ 1.62 (with z* ~ 1.97). At 95% c.l. we also obtain 
0.27 < fl m < 0.37 and 0.68 <h< 0.72, which seem to be 



in a better agreement with the current accepted values for 
both clustered matter density and Hubble parameters. 

V. CONCLUSIONS 

We have discussed a decaying vacuum scenario which is 
indistinguishable from the standard model with a genuine 
cosmological term in what concerns the general features 
of the predicted cosmic evolution. The early radiation 
phase in this A(i)CDM model is unaffected by the pro- 
cess of vacuum decay, as well as the physical phenomena 
taking place at early times (e.g., the primordial nucle- 
osynthesis). The following era is dominated by dust for 
a long time, and only recently the varying cosmological 
term has become important. 

Here, we have presented some quantitative results 
which clearly show that, even in the current stage of the 
Universe evolution, our decaying vacuum scenario is very 
similar to the standard one. We have also statistically 
tested the viability of the model by using the most re- 
cent SNe la observations, as given in Refs. @,0|. For the 
so-called gold sample, we have found 0.34 < f2 m < 0.44 
and 0.62 < H < 0.66 at 95.4% c.l., with the reduced 
X 2 — 1.15, which is similar to the one obtained in Ref. 
for a flat ACDM model. The SNLS data in turn provide 
0.27 < ft m < 0.37 and 0.68 < H < 0.72 (with X 2 ^ 1-0), 
which is in better agreement with the currently accepted 
estimates for both parameters. From the analysis pre- 
sented above, we also noted that a more precise determi- 
nation of the transition redshift zt from upcoming SNe 
la data may be able to distinguish this scenario from the 
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standard model since their predictions for this quantity 
are considerably different [see Eqs. i|2"T|) and 

We also emphasize that an important observational as- 
pect that deserves a careful investigation concerns the 
growth of density perturbations in the realm of this 
A(t)CDM model. In this regard, a preliminary analysis 
indicates that the evolution of the matter contrast Sp/ p 
shows no considerable difference relative to the standard 
ACDM case. A complete study on the formation of large- 
scale structures in this class of scenarios will appear in a 
forthcoming communication. 

Finally, it is worth observing that, from the theoretical 
viewpoint, it would be interesting to investigate possible 
relations between our approach and other quintessence 
models or modified gravitational theories recently dis- 
cussed in the literature. For this purpose, we note that 



some authors have already pointed out a mathemati- 
cal equivalence between dark energy and scalar-tensor or 
other forms of ideal fluid with inhomogeneous equations 
of state (see, e.g., 0). 
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